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1. Introduction 

The study of the physical models on the Sierpinski carpet has been the focus of 
several works in the last two decades. However, most of these papers are based on 
approximate renormalization group techniques [21 E] ? on perturbative expansions |3] or 
on numerical simulations E! • Thus, very few rigorous results are known for this self- 
similar structure. Indeed, the Sierpinski carpet is an infinitely ramified fractal, where 
it is not possible to apply the exact decimation techniques, which allow to solve many 
models, defined on finitely ramified fractals, such as the Sierpinski gasket [HIH]. 

A relevant theorem is given by an upper and a lower bound on the return probability 
of a simple random walk [9 . These bounds prove that the spectral dimension ^0] 
of the Sierpinski carpet is smaller than two. Hence, from the generalized Mermin- 
Wagner theorem ^T] , the continuous symmetry spin models do not present spontaneous 
magnetization. On the other hand, there are no general criteria relating the behaviour of 
the discrete symmetry models, in particular Ising, with a simple parameter, describing 
the topology of the space where the model is defined. The main known results regard 
the euclidean lattices, where spontaneous magnetization is present if the dimension is 
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Figure 1. The first 4 generations of the Sierpinski carpet graph. The sites of the 
graph (denoted with points) represent the spins and the links (denoted with hnes) 
represent the interactions. 

> 2, and the exactly decimable fractals, which never present spontaneous magnetization 

m 

However, the approximate calculations |2l HI and the numerical simulations [H] 
suggest that an Ising spin system on the Sierpinski carpet is spontaneously magnetized 
at finite temperature. Recently, this property has been proved in |il2j, where it is shown 
that it follows from the behaviour of percolation on the Sierpinski carpet. Here we give 
an alternative proof inspired by the classical Peierls argument ^H] for the Ising model 
on the two-dimensional lattice and, in particular, by the version given in ^3]. These 
theorems are the first exact results showing the existence of Ising transition in a low 
dimensional structure, i.e. in a structure with spectral dimension smaller than 2. 

We will consider the Sierpinski carpet graph illustrated in Figure 1 , even if the proof 
can be generalized to different infinitely ramified carpets. To each site of the graph we 
associate an Ising spin and to each link a ferromagnetic interaction. We call p~ the 
probability for the magnetization of the spin i to be smaller than one and we show that, 
at low enough temperature, p^ < e < 1/2 for each site i. The theorem is based on a 
low temperature loop expansion of this probability and on an estimate of the number 
of the self-avoiding on the sites (SOS) loops. 

The main difference from the two-dimensional case is that now the space, called 
dual graph, where the loops are defined, has unbounded connectivity, namely, the 
coordination (number of nearest neighbours) of a site is finite but not bounded. Hence 
the proof requires more complicated mathematical techniques. Notice that the standard 
Peierls argument results rather simple, as the dual of the two-dimensional lattice is again 
a two-dimensional lattice. In spite of the unbounded connectivity, here we will be able 
to estimate the number of SOS loops because of an important property of the dual 
graph, which will be called sparse high connectivity, i.e. sites with large coordination 
are far apart. In particular, when the coordination is greater than 4-3'", the sites are 
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distant at least 2"^+^ in the intrinsic metric of the dual graph. 

Peierls argument is valuable because of the insight it gives into why the Ising model 
is magnetized. Hence, we believe that this proof could be useful to find a general criterion 
for the existence of spontaneous magnetization on a generic discrete structure. Notice 
that Peierls argument has been already used to show the existence of Ising transition 
in an inhomogeneous systems, in particular in the case of diluted random graphs above 
the percolation tresholds We start introducing some definitions and notations. 



2. Definitions and notations 



Let us consider the generation of a Sierpinski carpet. On each site i we define an 
Ising spin Si = ±1 and we fix the spins on the boundary to the value +1 (see Figure 2). 
The energy of the spin system is given by: 

- ^ E ^i^i (1) 

where the sum runs over all the pairs (ij) of nearest neighbour sites, i.e. sites connected 
by a link. In the canonical ensemble the probability for the magnetization of the spin i 
to be smaller than 1 is: 



SiS-j 



{si,...sjv} 

In the denominator we sum over all possible spin configurations with +1 boundary 
conditions. In the numerator over all configurations with +1 boundary conditions and 
the spin i fixed to —1. We will prove the spontaneous magnetization of the system 
by showing the existence of a temperature (3 such that, in the thermodynamic limit, 
yp > (3 p- < e < 1/2. 

The theorem will be proved in the following steps. In Section III we introduce the 
dual space and the loop expansion of equation (j2)). With this expansion, we prove that 
p^ can be overestimated by the sum of exp(— over all SOS loops enclosing i (|/| is 
the length of the loop /). 

In section IV we introduce a tree representation of the SOS loops of length \l\ 
enclosing i and passing the site a. This representation allows us to overestimate the 
number of the loops as Zja YiaeF" ^a- 1" is a SOS loop. The product runs over the 
sites belonging to If excluding the site of highest coordination. Zja is the coordination 
number of the site representing, in the tree, the site of highest coordination in If. 

In section V we define the sparse high connectivity property and, using this 
property, we show that HaGZ"'" — ^^P(-^KI) where i^' is a suitable real constant. 

In section VI we prove that Zja < 10|/P°S2(3)_ particular we exploit the fact 

i 

that Zja is smaller than the number of the nearest neighbours, of the site of highest 

i 

coordination, which can be crossed by a loop of length |/| and surrounding i. 
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Figure 2. A spin configuration of the Ising spin model defined on the fourth generation 
of the Sierpinski carpet. The spins of the boundary are fixed to be +1. With the dotted 
lines we represent the interactions. 



In section VII we show that the number of sites, which can be crossed by a SOS 
loop of length \l\ and surrounding i, is smaller than 7r(8|Z|)^^°^2{3)^ This last estimate is 
proved by evaluating the maximum euclidean distance that can be covered in the dual 
lattice by a loop of |/| steps. Finally we will obtain 

oo 

P7 < E C|Z|3i°s2(3)g{K-/3)K| (3) 

|i|=0 

where C is a suitable positive constant. From equation Q we have that p~ < e < 1/2, 
for small enough temperature. 

Beside the sparse high coordination property, in the proof we strongly use the fact 
that the dual lattice is embedded in a two dimensional euclidean space and the fact that 
the plaquettes of the carpet are convex sets of this space. 

3. Loop expansion 

In analogy with Peierls argument, we write the Hamiltonian of the system in terms of 
a loop expansion. These loops belong to the dual graph of the discrete structure where 
the model is defined. In particular, the dual graph of the Sierpinski carpet is obtained 
by introducing a site a for each plaquette of the carpet. Then, two sites are connected 
if and only if it exists a link of the carpet separating the plaquettes relevant to the sites 
we are considering. Hence, there is a one to one correspondence between the links of the 
carpet and the links of the dual. We will call the set of the sites of the dual graph 
and Ed the set of the links. Notice that the dual graph has unbounded connectivity; 
hence, it must be studied with a certain caution. The length (number of links) of the 
shortest walk connecting two sites G Vd will be called the chemical distance between 
them. In Figure 3 we illustrate the Sierpinski carpet and its dual graph. 
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Figure 3. The dashed Une denotes the Sierpinski carpet and the continuous hne the 
dual graph. The site i belongs to the carpet, and a to the dual graph. 



Let us introduce a one-to-one correspondence between the spin configurations and 
the self-avoiding loops of the dual graph (self-avoiding on the links). A spin configuration 
with +1 boundary conditions is identified by the links of the carpet which separate the 
spins with -|-1 magnetization and the spins with —1 magnetization. Since, by definition, 
each link of the dual graph corresponds exactly to a link of the carpet, we obtain a 
relation between the spin configuration and a subset L C E^. It is easy to show that 
the number of links in L adjacent to the same site a is even. Therefore, L can always 
be obtained as a superposition of self-avoiding (on the links) loops. If in L there are 
more than two links adjacent to the same site, the decomposition of L into loops is not 
unique. 

On the other hand, let us call L any subset C Ed such that for any a E Vd the 
number of links in L adjacent to a is even, i.e. L can be decomposed into loops (we 
will call El the set of all possible L). In a unique way we can associate to L a spin 
configuration with +1 boundary conditions, by fixing to +1 the spins enclosed by an 
even number of loops and to —1 the spins enclosed by an odd number of loops. In 
Figure 4 we illustrate a spin configuration and its corresponding L C Ed. 

The energy of the spin configuration can be written in terms of the corresponding 
set of links L. Indeed 

^ = E = -l^^l + \ E(^^ - ^.)' (4) 

where \Ed\ is the total number of links of the dual graph. A link gives a contribution 
1 to the Hamiltonian (j3)), if it connects opposite spins, and zero, if it connects parallel 
spins. Dropping the irrelevant additive constant \Ed\, we obtain H = \L\, where is 
the number of links in the set L G El corresponding to the spin configuration of energy 
H. If we consider a decomposition of L into loops / we obtain 

H=\L\=J2\l\ (5) 
I 
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where |/| is the length of the loop /. The probability p^ is 



E 



e 



-/3\L-\ 



Pi = ^ -m 



E 



e 



In the denominator we sum over all possible sets L G El and in the numerator over all 
subsets G such that an odd number of loops encloses the site i. E~[ denotes the 
set of all G El, such that an odd number of loops encloses i] while denotes the 
set of all L'^ G E^, such that an even number of loops encloses i. Hence, E^ = E^ U Ej^. 

We will call /j the shortest SOS loop belonging to L~ and surrounding i (see Figure 
5). Each set L~ can be considered as the union of U and a suitable set of E^. Therefore, 
equation (jH)) can be written as 



E e-^'" '+ E 



denotes the sum over all possible SOS loops k enclosing the site i. 
sum over all G E^ such that fl /j = and such that, in k U L"*", is the shortest 
SOS loop surrounding i. Dropping the restriction on the sum we obtain an 

overestimate of p~ 

m \^ g-/3|i+i 
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ni{k) represents the number of SOS loops of length k enclosing i. Eventually we can 
write 

oo 

pr<EE<(fcK^' (9) 

where n°'{k) is the number of SOS loops of length k, enclosing i and passing the site a 
of the dual graph. 

Now we take the thermodynamic limit. The theorem will be proved by showing on 
n°'{k) an exponential bound independent from i. Up to now we followed step by step the 
Peierls argument for the Ising model on two-dimensional lattices jT3] . However, here the 
dual graph is more intricate than in the two-dimensional case, since it has unbounded 
coordination. Hence, the proof of the bound on nf{k) requires more work. 

4. Tree representation of the loops 

In order to count the SOS loops of length k, surrounding i and starting from a, we 
represent them into a tree. By induction on we prove that any set of N SOS loops of 
length k, starting from a', can be represented in a tree satisfying the following properties 
(TP). 

Each site of the tree a is associated to a site of the dual graph a (a). Nearest 
neighbour sites in the tree correspond to nearest neighbour sites in the dual graph. 
The correspondence is not one-to-one, however, if b and c are nearest neighbour of the 
same site of the tree, then they represent different sites of the dual, i.e. a{b) ^ a{c). 
Consequently, calling Za the coordination number of a site of the tree and the 
coordination of a site of the dual, we have Za < Za(a)- One site of the tree is called 
the origin o. We have a{o) = a' and the maximum distance of any sites of the tree from 
is k {k will be called the length of the tree). There are exactly sites at a distance 
of k steps from a and to each of them we associate, in a one-to-one correspondence, one 
of the loops. These sites will be called ending points e. a{e) = a' and they are the only 
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ones in the tree with coordination 1. Therefore, the number of loops represented in the 
tree can be evaluated by counting the ending points. Finally, the sites crossed in the 
tree by the walk of k steps, connecting the origin to an ending point e, correspond, in a 
one-to-one relation, to the sites of the dual graph belonging to the loop associated to e. 

The set made of any single loop can be represented by a linear chain of k + 1 sites 
and k hnks. The first site of the chain is the origin and it corresponds to a', the second 
corresponds to the second site of the loop and so on, the last site corresponds again to 
a' and it is the ending point of the tree. 

Let us show that, if there is a tree representation for any set of — 1 loops starting 
from a', we can add to any of these sets a new loop and produce a tree representing this 
new set. We will proceed by adding a new chain to the tree, relevant to the A^ — 1 loops. 
The sites of this tree correspond, in both representations, to the same sites of the dual 
graph and they represent the first — 1 loops. The starting point a' of the new loop 
is naturally associated to the origin a of the tree. Then, if one of the nearest neighbour 
of corresponds in the dual to the second site of the loop, we move to this site of the 
tree and check if one of its neighbours correspond to the third site of the loop and so 
on. Until we find that, at the step q, no neighbour of the site represents the site reached 
in the q + 1 step of the loop. Since we are considering SOS loops (we can not move 
backward), this site is distant q from the origin of the tree. Now we add to this site a 
bifurcation and a new chain of length k — q. The first site of this new chain corresponds 
to the site reached in the g + 1 step of the loop, the second to the site reached in the g + 2 
step, the last corresponds again to a' . So, when we add a loop to the set, we add to the 
tree one ending point naturally associated with the new loop. The new neighbour sites 
added to the tree are neighbour also in the dual. Furthermore, we never add to any site 
of the tree a a new neighbour 6, such that a{h) is already represented in the tree by a 
neighbour of a. This way, we produce a tree satisfying TP properties and representing 
a set of N loops passing a' . 

In particular, it is possible to construct the tree representing the set of all SOS loops 
of length k surrounding i and starting from a. The number of these loops is equal to the 
number of the ending points in the tree. In Figure 6 is illustrated a tree corresponding 
to the all the loops of length 7 passing a = and enclosing i. 

Let us consider a tree of length k and satisfying TP. We call e the ending point 
maximizing the product Pk of the coordination numbers of the sites crossed by the k 
steps walk connecting the origin to e. We prove by induction on k that the number of 
the ending points is smaller than P^. 

For k = 1 the tree is given by the origin and its Zg neighbours which are the ending 
points. Hence both Pi and Ti are equal to Zo- Let us now suppose that Tj^-i < Pk-i for 
any tree of length k—1. We need to show that, given a tree of length k, < P^. Cutting 
the origin o from the tree of length k, the graph is disconnected into Zg subgraphs. Each 
of these subgraphs is a tree of length k—1, satisfying TP, whose origin is one of the 
neighbours of o. We will denote these trees with the index h and we will call T^_i the 
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number of the ending points of the tree h. is equal to the sum over h of T^_i. Hence: 

n = f: Tt, < f: Pt, < m^{Pt,) (10) 
h=l h=l 

In the first inequahty we use the fact that Tj^_^ and Pk_i are, respectively, the number 
of the ending points and the maximum of the products of the coordination numbers in 
a tree of length k — 1. Since the origin belongs to any of the walks connecting a to an 
ending point, Pk can be expressed as the product of Zo and the maximum value of -P^_i. 
Therefore, from one immediately obtains the proof. 

In the tree representing all the SOS loops of length \l\ enclosing i and passing a, we 
will call wf the walk connecting the origin to e and we will call If the loop associated 
to e. We have: 

Km) < n = % n < %exp ( ^ \og{zM (n) 

Zja is the coordination number of the tree site m corresponding to the site u of the loop 

_ i 

If with the highest coordination in the dual graph (if the sites of highest coordination 
are more than one, we arbitrarily chose one of them), w'" and are obtained by 
subtracting from the walk wf and the loop If, respectively, m and /i. In the last 
inequality we use the one-to-one correspondence between w'" and P° and fact that the 
coordination number Za of the sites of the tree is smaller or equal than the coordination 
number z.y(^a) of the relevant sites of the dual graph. 
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5. Sparse high connectivity 

In the dual graph the possible values of the coordination numbers are 4 ■ 3"^ with 
m = 0, 1, . . .. Hence we have 

<(|/|) < %exp (^J2Npt{m){\og{4)+m\og{3))^ (12) 

Np'^ (m) is the number of sites in with coordination number in the dual graph equal 
to 4 ■ 3"^. 

The minimum chemical distance d{m) between two sites of the dual graph with 
coordination greater or equal than 4 ■ 3*" is: 

d{m) = l \, '1"^ = ^ (13) 
V J I 2™+! if m > ^ ^ 

Equation (fT^ represents the fundamental topological property of the dual graph of the 
Sierpinski carpet which will allow us to prove the exponential bound on ni{k), even in 
absence of bounded coordination. We will call this property sparse high connectivity, 
since it implies that, in the topology of the dual graph generated by the chemical 
distance, sites with high coordination are far apart. 

For a SOS loop /, the sites with coordination > 4 ■ 3™ are less than 1 + \l\/d{m). 
Since in T'" we have excluded the point of highest coordination, we can overestimate the 
number of sites of coordination 4 ■ 3"^ as: 

< (14) 

From (fT^ we obtain that, if m > log2(/) — 1 = mM, then Npo'{m) = 0. Therefore, tum 
represents the maximum value of m on the sites of P". 
Thus we have: 

\l\ log(4) + / 2-™-'(log(4) + mlog(3)) j 

< % exp [|/| log(4) + / £ 2— i(log(4) + m log(3))] 

< Zfo. exp {K\l\) (15) 
where is a suitable positive, real constant. 



6. The site of highest coordination 

In this section we will show that Zj^ < 10/^°S2{3)_ if < io/i°g2{3) jg the site in the 
dual graph with highest coordination), the inequality is trivially satisfied. Otherwise, 
we will use the fact that Zja is the coordination of a site of a tree representing SOS 

i 

loops surrounding i. Two sites of the tree are neighbour only if the corresponding sites 
of the dual are neighbours and it exists a SOS loop passing through them and enclosing 
i. Therefore, the bound on Zja is proved by showing that the Z j. < 101^"^^^^^; where 



Ising model on the Sierpinski carpet fractal 



11 




Figure 7. The plaquette of the site is given by the region denoted by the dashed 
Une. The sites within this region (excluded fj,) are the sites belonging to the plaquette. 
The distance along the boundary between the sites of the plaquette i' and p is the 
length of the line plotted with the heavy line. In this example fp.^ — 6. With the 
dotted lines we denoted a walk in the dual graph between v and p of length 18. 



Z^j. represents the number of nearest neighbours of /i, which can be crossed by a SOS 
walk of length |/|, passing /x and enclosing i. 

Let z/ and p be nearest neighbour of p. We define the distance along the boundary 
f^^p as in Figure 7. In Figure 7 the plaquette of the site p is also defined. 

Let us consider a SOS loop l'^'^ of length |/| passing p and its nearest neighbours p 
and u. We will prove that: 

Furthermore, if l'^''' encloses the site of the carpet i, from the planarity of the dual graph 
we have the following property: no loop of length \l\ passing p and enclosing i can pass 
through any site r of the plaquette of p, such that r is located with respect to p in the 
opposite direction of z/ at a distance f^^p > 4 ■ /'°S2(3)_ 

An analogous property holds for the sites of the plaquette at a distance greater 
than 4 ■ /'°S2(3) from u in the direction opposite to p. Therefore, since f^^p < 2|/p°^2(3)^ 
the number Z^j. of sites belonging to the plaquette of p and crossed by a loop passing 
p and surrounding i is smaller than 10|/|^°^2(3)_ This way we obtained the bound on Zja. 

i 

Let us prove ()16|) . The dual graph is naturally embedded in a two dimensional 
euclidean space. Since the plaquettes are convex subsets of this space, the length of the 
border f,y^p turns out to be smaller than the length, in the euclidean space, of any walk 
from u to p which do not pass through the plaquette itself. In particular, it is smaller 
than the length of l'"'''. l'"''' is the walk obtained subtracting p from the loop l'^'''. Now, 
a step in the dual graph, crossing a site with coordination Zp, covers, in the euchdean 
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space, a distance smaller than Zp/2. We have: 

1 ruM log2(l'l)-l /Q\m 

ru,p < - E -^P < 2 E Ni>.Am)3'^ < 21 + I E ( o ) < 2/'°'^^'^ (17) 

Ni'i^,p{m) is the number of sites in 1"^'P with coordination 4-3™. In equation p7|) Niv.p{m) 
has been estimated using inequahty ((T3j). 

Now we need to prove that no that no loop of length |/| passing r and /i can enclose 
the site i. We suppose, ad absurdum, the existence of this loop and we call it V''^ (a is 
the other site of the loop belonging to the plaquette of /i) . Since r and a are connected 
by a loop of length from we have f^^a < 2|/|'°S2(3). Hence fp^„ > 2\l\^°^^^^\ If also 
r''^ surrounded i, for the planarity of the dual graph, there should exist two sites 7 and 
6 belonging both to Z"^'"" and to Z'^''' (see Figure 8 for an illustration). 

We denote with |7,5|;t,<t the length of the part of the walk within the sites 7 
and 6. We have 

|p, 6\i^.p + \6, 7|i,.,p + I7, u\i.,p = \l\-2 

\a, (5|ir,a + \S, 7|ir,<T + I7, T\lT,a = \l\ - 2 (18) 

and then: 

|p, S\i.,p + |7, u\i..p + \a, S\ir,. + |7, r|r- < 2{\l\ - 2) (19) 
On the other hand, since fj,,- > |/|'°S2(3) Q^Yid fp^^^ > \l\^°^'2^^\ from ()16|) we get 

|p, 5\li^,P + |(T, 5\lr,a > \l\ -2 

|z/,7|i..p + |r,7|zr.. > |/| - 2 (20) 

and then 

|p, 6\i^,p + \a, 6\ir,. + \v, 7|z.,p + |r, 7|r,- > 2(|/| - 2) (21) 

Equation (j2H) can not hold together with equation (jl9j) . Therefore, assuming the 
existence of the site r, we obtained the contradiction we wanted to prove. Hence, 
we have concluded the proof of the bound on Zja. 



Ising model on the Sierpinski carpet fractal 



13 



7. Euclidean length of a loop 

With the bound on Zfa proved in the previous section, we obtain that probability for 

i 

the site i to have negative magnetization can be estimated as 

oo 

p- < E Ari(|/|)10|/|'°S2(3)e(^-/3)l^l (22) 
\i\=o 

iVj(|/|) is the number of points in the dual lattice which can be crossed by a SOS of 
length \l\ which encloses i. 

An estimate of A^i(|/|) can be obtained by using again the fact that the dual graph 
is embeddable in a two-dimensional euclidean space. In particular, in a walk in the 
dual graph, a step crossing the site a covers in the euclidean space a distance smaller 
or equal than z'^ = fa-^a+, where a~ and a"*" are respectively the sites preceding and 
following a within the walk. Let us consider /j SOS loop in the dual graph of length |/| 
and enclosing the site i of the carpet. We call d{li) the distance covered by this loop in 
the euclidean space, we get 

d{h) < E < < + E ^« (23) 

In the sum we separate the contribution of the site with highest coordination number 
and we used the simple property that z'^ < Za- Zj represents the distance along the 
boundary between the sites preceding and following the site of highest coordination 
number in a walk of length |/|. In the previous section we proved that < 2/'°S2(3)_ 
Furthermore the sum over all sites belonging to /• has been already evaluated in (|T7j). 
Therefore we have that 

d{li) < 8|/|^°S2(3) (24) 

The number of sites of the dual graph, distant from i in the euclidean space less than 
d{k) it is smaller than nd^ik). Therefore we obtain that Ni{l) < 7r{8\l\^"^^^^^f. Hence 

oo 

p- < E ^25) 

\l\=0 

where C is a suitable real constant. For small enough values (3 the sum in ()25|1 is 
smaller than e < 1/2. Thus we obtain the uniform bound on p~ proving the existence 
of spontaneous magnetization at finite temperature. 
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